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^ Accurate  two-stream  representation  of  the  angular  distribution  of  radiation  exchange  in  a plane- 
parallel  atmosphere  requires  the  use  of  a diffusivity  factor  at  least  1 3%  larger  than  the  normally 
chosen  values  of  1.66  or .732},  since  the  flux  divergence,  rather  than  the  transmitted  flux, 
determines  the  heating  ocfeoTing  rates.  Diffusivity  factors  are  calculated  to  yield  exact  results 
at  various  optical  depth/for  Doppler  and  Lorentz  line  profiles  as  well  as  for  the  grey-absorption 
case. 
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TWO-STREAM  PARAMETERIZATION  OF  THE 
FLUX  DIVERGENCE  IN  A PLANE-PARALLEL  ATMOSPHERE 

1.  INTRODUCTION 

To  speed  up  the  solution  of  time-consuming  radiative  transfer  problems  in  plane-parallel 
media  the  technique  of  the  two-stream  approximation  is  often  adopted,  particularly  for  plane- 

■ ' , I \ 

tary  atmospheres.0  2-31  The  reciprocal  of  the  cosine  of  the  angle  that  represents  the  mean  obli- 
quity of  the  inward  and  outward  streams  is  known  as  the  "diffusivity  factor,”  dr.  The  two  most 
common  choices  for  its  value  are  1.66°  21  and  >/3  (—  1 .732) .<4>  The  value  of  1.66'5'  follows 
from  the  fact  that  a monochromatic  isotropic  radiation  field  is  attenuated  by  one-half  at  a 
monochromatic  normal  optical  depth  of  0.42,  that  is,  0.5  — exp  (—  1.66  x 0.42).  The  n — 2 
Gaussian  quadrature  integral  approximation,4)  yields  a value  of  >/3,  which  differs  from  1.66  by 
only  4.3  per  cent. 

Two  important  atmospheric  situations  render  these  values  inaccurate.  First,  in  many  cases 
of  physical  interest,  the  radiation  field  is  highly  anisotropic.  Second,  net  exchange  of  radiative 
energy  (as  opposed  to  simple  transfer  of  radiation  from  one  layer  to  another)  can  occur 
between  layers  separated  by  optical  depths  much  smaller  or  much  larger  than  0.42.  For  exam- 
ple, radiation  exchange  in  the  earth’s  atmosphere  by  the  15/a  CO2  bands  generally  occurs 
between  layers  separated  by  large  optical  depths  since  the  atmosphere  is  isothermal  on  a scale  of 
only  0.42  optical  depths  (less  than  30  m at  an  altitude  of  35  km).  It  should  be  emphasized  that 
the  flux  divergence,  rather  than  the  transmitted  flux,  determines  the  net  exchange. 
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2.  BASIS  FOR  CHOOSING  THE  DIFFUSIVITY  FACTOR 

Consider  two  differentially  thin  atmospheric  layers  at  z'  and  z".  Monochromatic  radiation 
emitted  from  the  layer  at  z'  toward  the  layer  at  z"  has  an  angle-averaged  probability  P,  of 
reaching  z"  given  by 


So  A'W  «p  ("  -1)  d» 

f dn 


(1) 


In  equation  (1),  n is  the  cosine  of  the  angle  of  obliquity,  jHn)  is  the  emission  coefficient  at  z', 
and  rI  >Z"  is  the  monochromatic  normal  optical  depth  between  z'  and  z".  If  jt  is  isotropic,  i.e., 
not  a function  of  m.  we  find 


Pr  - Jo  CXP 


T . 

2 .2 


dfi  - Ej,irz  y) 


(2) 


where  £2  is  the  familiar  second  exponential  integral. 


Clearly,  the  amount  of  radiation  that  originates  at  z'  and  is  absorbed  in  a layer  dz  at  z"  is 
proportional  to  dPe/dz  at  z".  Let  Pe(r)  be  the  monodirectional  escape  probability  in  the  perpen- 
dicular direction  (/a  — ± 1).  If  a diffusivity  factor  df(r)  can  be  found  which  correctly 
represents  P,(.r)  by  Pf(dfr),  the  angular  dependence  of  the  radiative  transfer  will  be  accurately 
parameterized  by  dr.  The  fundamental  importance  of  Pe(r)  has  been  pointed  out  by  Dickin- 


son 


(6) 


The  quantity  P,(r)  is  also  equal  to  the  factor  by  which  the  flux  fis  diminished  across  the 
optical  path  r,  where 

F “ 2ir  S_t  d/i 
2 
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and  /(/*)  is  the  specific  intensity.  The  angular  distribution  of  /,  i.e.  its  m dependence, 
significantly  affects  the  value,  of  P,  across  a given  path.  If  I (pi)  “ /,  (isotropic)  across  a hemi- 
sphere, then,  across  path  tq. 


P&J  - F(t„)/F(t  - 0) 


■ ..  : >t  i i it 


X w (-  ‘tt) 


— I p.dpL/do/2) 


■ 2Ej(to),  . (3) 

Such  an  isotropic  flux  would  be  produced  by  a Lambert  surface.  However,  the  flux  emitted 

• . j ’ f * » 

from  an  optically  thin  layer  of  optical  depth  A in  the  atmosphere  is  characterized  by  a highly 
aniSQtropic - I (pi)-  In  fact,  /(ft)  - SA(fi)  - SA/ft  where  S is  the  source  function  j/k.  The 
corresponding  fluxdirninutipn  P,(*a)  is  equal  to  , 

F(r0)/F(0)  '*  /0'('SA/m)  pcdpl J0'(SA/mWm 

t—  E2(*<>).  'j  • .» . W) 

precisely  the  result  of  Equation  (2). 'the  physical  content  of  equations  (2)-(4)  may  be  summar- 
ized as  follows.  The  escape  probability  (flux  transmission)  across  a layer  equals  2E3  if  the  radi- 
ation incident  upon  the  layer  is  isotrbpic.  Heating  or  cooling  due  to  exchange  between  layers  is 
determined'  by  the  transmission  of  radiation  among  differentially  thin  layers  whose  escape  pro- 
bability across  a finite  layer  equals  E2 . The  fact  that  the  integrated  radiation  field  from  an  opti- 
cally thick  finite  portion  of  the  atmosphere  may  be  nearly  isotropic  is  not  relevant  since  it  is  the 
derivative  of  the  2 £3  transmission  function  («  E2 ) for  isotropic  radiation  which  determines  the 
exchange.  Thus  diffusivity  factors  based  upon  E2  rather  than  2 £j  should  be  used  whenever  the 
calculation  of  exchange  phenomena  such  as  heating  or  cooling  is  the  objective. 

If  radiation  is  transferred  primarily  within  lines  the  most  useful  quantity  to  consider  is  the 
line-profile  averaged,  angle-averaged  escape  probability 
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Pr^o,,-  f*~  *(x)  Pr(rx)dx  (5) 

In  equation  (3)  jt  is  the  frequency  variable,  4>(x)  is  the  line  emission  profile,  assumed 
proportional  to  the  absorption  profile,  and  normalized  such  that 

J m ^(x)  dx  — 1 (6) 

The  subscripts  d or  / correspond  to  the  P,  for  a Doppler  or  Lorentz  line  profile,  respectively. 

For  a Doppler  profile  <f>(x)  - — for  a Lorentz  profile  *(x)  - 1/»(1  + x2). 

3.  NUMERICAL  RESULTS 

The  integrals  of  Eqs.  (2)  (for  monochromatic  radiation)  and  (5)  (for  both  Doppler  and 
Lorentz  profiles)  were  computed  numerically  and  compared  with  the  monodirectional  escape 
probabilities  to  find  the  difTusivity  factor  required  for  exact  agreement  with  Pr.  In  Fig.  1 the 
line  profile-averaged  escape  probabilities  are  plotted  for  Doppler  and  Lorentz  lines  for  optical 
depths  from  10-1  to  I0J.  Figure  2 displays  the  required  diffusivity  factors  up  to  7 optical  depths 
at  line  center,  the  most  interesting  region  of  variation.  Also  shown  in  Fig.  2 is  the  diffusivity 
factor  for  monochromatic  radiation  (frequency-independent  absorption).  For  7 < r < 103, 
df  — 1.95  for  Doppler  lines,  while  the  df  increases  from  2.22  at  r — 7 to  2.24  at  r - 103  for 
Lorentz  lines.  For  grey  absorption  (monochromatic  radiation)  df  slowly  approaches  1 as 

r — oo. 

The  limits  of  df  as  r -•  may  be  straightforwardly  derived  from  equation  (5).  Integrals 
similar  to  equations  (2.21)  and  (2.27)  of  Holstein*73  are  obtained  and  easily  solved  by  his  tech- 
nique as  r — * °*.  The  Doppler  diffusivity  factor  slowly  approaches  2 as  r - » while  for 

Lorentz  lines  the  limit  is  df  - jyj  - 2.25.  For  r - 0,  a diffusivity  factor  based  on  flux 
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transmission  ■ is  2.0(^,  whereas  a diffusivity  factor  appropriate  for  flux  divergence  approaches 
+oo.  The  latter  indicates  that  only  that  portion  of  the  radiation  field  which  traverses  a long  opt- 
ical path  (i.e.,  large  oblique  angles)  can  cause  heating  or  cooling  as  t — • 0. 


4.  SUMMARY  AND  CONCLUSIONS 


If  the  primary  objective  of  a radiative  transfer  calculation  is  to  obtain  radiative  exchange 
rates  rather  than  flux  transmission,  a diffusivity  factor  must  be  chosen  to  match  the  flux  diver- 
gence rather  than  the  flux.  It  has  been  shown  that  these  factors  for  an  optically  thick  atmo- 
sphere should  range  from  1.9S  for  a pure  Doppler  profile  to  2.2S  for  a pure  Lorentz  profile. 
These  values  are  significantly  higher  than  the  usual  value  of  1.66  for  heating  calculations. <2> 
The  need  for  larger  diffusivity  factors  has  been  confirmed  by  the  multiple  scattering  continuum 
heating  calculations  of  Lacis  and  Hansen, (3)  who  found  that  df  — 1.9  provides  the  best  fit  to 
their  detailed  numerical  work.  This  value  is  exactly  the  diffusivity  factor  required  for  optical 
depth  unity  and  grey  absorption  (Fig.  2),  where  the  bulk  of  continuum  transfer  effects  occur. 
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Figure  I — Line-profile-averaged  escape  probabilities  arc  logarithmically  plotted  against  optical  depth.  For  pure  Doppler 
and  Lorentz  profiles.  Values  are  given  For  the  connecting  path  perpendicular  to  the  planes,  and  For  the  angle-averaged 
escape  probability. 
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Figure  2 — The  required  diffusivity  factor  is  given  for  line-profile-averaged 
escape  probabilities  for  Doppler  and  Lorenlz  lines,  and  monochromatic  radia- 
tion, as  a function  of  optical  depth. 


